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Abstract. Dense wind of a massive star can be partially captured by a neutron star (NS) inside a compact binary system. 
Depending on the parameters of NS and the wind, the matter can penetrate the inner NS magnetosphere. At some distance from 
the NS a very turbulent and magnetized transition region is formed due to the balance between the magnetic pressure and the 
pressure inserted by accreting matter. This region provides good conditions for acceleration of particles to relativistic energies. 
The matter at the transition region can farther accrete onto the NS surface (the accretor phase) or is expelled from the NS 
vicinity (the propeller phase). We consider the consequences of acceleration of electrons at the transition region concentrating 
on the situation in which at least part of the matter falls onto the NS surface. This matter creates a hot spot on the NS surface 
which emits thermal radiation. Relativistic electrons lose energy on the synchrotron process and the inverse Compton (IC) 
scattering of this thermal radiation. We calculate the synchrotron spectra (from X-rays to soft y-rays) and IC spectra (above 
a few tens MeV) expected in such a scenario. It is argued that a population of recently discovered massive binaries by the 
INTEGRAL observatory, which contain neutron stars hidden inside dense stellar winds of massive stars, can be detectable by 
the recently launched Fermi LAT telescope at GeV energy range. As an example, we predict the expected y-ray flux from 
recently discovered source IGR 119140+0951. 
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1. Introduction 

Massive binary systems, containing a compact object (neutron 
star or black hole), have been suspected as sites of particle ac- 
celeration up to at least TeV energies since early 80'ties. But 
only recently, three such type of objects have been confirmed 
as VHE -y-ray sources by the Cherenkov telescopes. It is clear 
that at least in one case the y-ray emission is due to accelera- 
tion of particles as a result of interaction of the energetic pulsar 
wind with the wind of a massive star. Other TeV y-ray binaries 
either operate under this same general scenario or the accel- 
eration of particles occurs in the jet launched by the compact 
object. In the case of these TeV -y-ray binaries, the neutron star 
(pulsar) should rotate fast enough in order to produce energetic 
pulsar wind. 

In this paper we consider the case when the neutron star 
is not so energetic and the matter from the stellar wind can 
penetrate into the inner pulsar magnetosphere. At certain dis- 
tance from the NS surface the magnetic pressure balances 
the pressure of accreting matter creating very turbulent tran- 
sition region. Then, the NS is enshrouded in a dense co- 
coon and resembles objects called as the hidden pulsars (see 
e.g. Tavani 1991, Tavani & Brookshow 1993) or cauldrons 
(Begelman & Rees 1984, Treves et al. 1993). Note that y- 
ray emission has been predicted from such hidden pulsars 



(Tavani 1993). We consider the situation in which electrons are 
accelerated in turbulent, strongly magnetized, transition region 
in the inner pulsar magnetosphere. In principle, electrons can 
even reach ~ TeV energies but, as we show below, due to large 
synchrotron energy losses only y-rays in the Fermi LAT tele- 
scope energy range may be effectively produced. 

It looks that the INTEGRAL observatory has recently dis- 
covered a class of objects within the sample of compact high 
mass X-ray binaries (see e.g. Chaty et al 2007, Rodriquez & 
Bodaghee 2008) which might operate under the mechanism 
discussed by us. These newly discovered massive binaries are 
compact, with the orbital periods a few to several days. Some 
of them contain relatively slowly rotating neutron stars which 
may allow the matter to penetrate the inner NS magnetosphere. 
According to the classification scheme of X-ray binaries (e.g. 
Lipunov 1992) such binaries with slowly rotating neutron stars 
(NS) belong to the propeller and accretor class. In these objects, 
the matter of the dense stellar wind accrete onto the strongly 
magnetized NS, interacting with the rigidly rotating inner mag- 
netosphere. In the case of accretor most of the matter reach the 
NS surface. In the case of propeller most of the matter is ex- 
pelled from the NS vicinity due to the centrifugal force. 

It is proposed that such general scenario can be applied to 
a type of obscured compact objects inside massive binaries re- 
cently observed by the INTEGRAL observatory. We perform 
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detailed calculations of the y-rays spectra including in some 
cases also the inverse Compton (IC) e* pair cascade initiated by 
relativistic electrons in the radiation field produced on the neu- 
tron star surface. We also take into account the synchrotron en- 
ergy losses of primary electrons and cascade e* pairs, in order 
to obtain simultaneous synchrotron X-ray spectra. The acceler- 
ation of hadrons in such a scenario and possible production of 
high energy y-rays, neutrinos and neutrons will be discussed in 
another paper. 

2. Description of the model 

We consider a compact binary system containing rotating neu- 
tron star (NS) and a massive companion of the 0,B type star. It 
is assumed that a mass from the stellar wind is effectively cap- 
tured by a strong gravitational potential of the NS. Such slowly 
rotating neutron stars appear at a certain stage of the evolution 
of the binary system due to the loss of angular momentum ei- 
ther by the pulsar mechanism or as a result of the torque exhib- 
ited by the matter accreted from the stellar wind. Depending 
on the rotational period and surface magnetic field of NS, 
the accretion process onto NS can occur in different phases. 
According to the classification scheme of Lipunov (1992), the 
accretion process can occur in the phase of accretor (for rel- 
atively slower rotators) or in the phase of propeller. For very 
energetic, strongly magnetized and short period neutron stars, 
the accretion process does not occur at all since the matter can 
not penetrate below the light cylinder radius. In the case of NS 
in the accretor and propeller phases, the matter from the stel- 
lar wind penetrate below the light cylinder radius of the rotat- 
ing NS magnetosphere. This matter extracts rotational energy 
from the NS as a result of the interaction of a free falling matter 
with the rigidly rotating inner NS magnetosphere. In this paper, 
we consider such direct collision of the matter from the stellar 
wind with the NS magnetosphere (see general scenario shown 
in Fig. [TJ. As a result of such interaction, a very turbulent and 
magnetized transition region is formed. In the case of accretor, 
most of the matter falls onto the NS surface creating small hot 
region on the NS surface. In the case of propeller, most of the 
matter is expelled from the vicinity of the NS. It is expected 
that the matter might be able to accrete also in this stage but 
in a rather non-stationary form since due to gradual accumula- 
tion of the matter close to the transition region the pressure of 
the matter can overcome the pressure from the rotating mag- 
netosphere. It is supposed that some amount of the matter can 
accrete but most of it is expelled from the binary system (form- 
ing a larger scale jets ?) for some specific transition parameters 
describing the accretor/propeller scenario. 

Note that in some massive binaries (e.g. Her X-l), the ac- 
cretion of matter can occur also through the Roche lobe over- 
flow. Then, the accretion disk around the NS is at first formed 
and the matter from the disk interacts directly with the NS mag- 
netosphere. We do not discuss this possibility in this paper but 
rather concentrate on a simpler quasi-spherical accretion onto 
magnetized NS of the matter from the dense wind of the mas- 
sive star. The accretion through the disk can be easily taken into 
our considerations in our scenario by introducing a geometri- 
cal factor which describes the part of the sphere in which the 




Fig. 1. Schematic representation of the model discussed in the 
paper. The wind from a massive star penetrate the inner mag- 
netosphere of a neutron star. As a result of the interaction of the 
wind with the dipole magnetic field of the NS a strongly mag- 
netized turbulent region is created (thick solid curve) in which 
the matter gains energy being accelerated to the velocity of the 
rotating magnetic field lines. Electrons (e) are accelerated in 
this turbulent transition region and cool on the synchrotron pro- 
cess and on the comptonization of thermal radiation (marked by 
X), emitted from the stellar surface around the magnetic pole 
region. As a result, the first generation of IC y-rays is produced. 
Thermal radiation appears due to the gravitational energy re- 
lease of accreting matter onto the NS surface. If the optical 
depths for y-rays are large enough, they initiate the IC e* pair 
cascade in this same thermal radiation. 

accretion occurs. Therefore, the accretion through the disk will 
correspond to one of our discussed cases with effectively larger 
accretion rate occurring only inside a limited region determined 
by the thickness of the accretion disk. 

Let us at first consider the case of the NS in the phase of 
pure accretor. The accretion rate of the matter onto the neutron 
star(M acc = 10 I6 M 16 gs-') can be estimated from the observed 
(in some cases) thermal X-ray emission (Lx = 10 36 L36 erg s _1 ). 
These two values can by related for the known radius and the 
mass of NS (we assume R^s ~ 10 6 cm and Mns = 1-4M Q ), 

M acc » 5 x 10 15 L 36 gs" 1 . (1) 

This matter arrives onto the neutron star surface along the mag- 
netic field lines. The distance at which the magnetic field starts 
to dominate the dynamics of the matter (the Alfven radius) can 
be estimated by comparing the magnetic field energy density 
with the kinetic energy density of the matter, 

Z?i/8;r = pv t 2 /2, (2) 

where Z?a is the magnetic field in the inner neutron star mag- 
netosphere^ = M- iC cl(nR 2 A V{) is the density of accreting mat- 
ter, Vf = (2GMns/Ra)^ 2 is the free fall velocity of accreting 
matter, Ra is the Alfven radius, and G is the gravitational con- 
stant. The medium in the transition region is very turbulent and 
strongly magnetized providing good conditions for the acceler- 
ation of particles to high energies. Let us estimate the location 
of this region from the surface of the neutron star by applying 
Eq. [2] and assuming that the magnetic field in the neutron star 
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magnetosphere is of the dipole type, i.e. Bp, = Pns(Pns/Pa) 3 , 

fi A =4x 10 % B^M~1 P cm, (3) 

where the magnetic field at the neutron star surface is Z?ns = 
10 12 Z?i2 (see also Baan & Treves 1973). Based on the known 
vale of Pa, we can estimate the magnetic field strength at the 
transition region, 

B A = 1.6 X 10 4 M^ 7 B n 5/7 G. (4) 

Considered here accretor scenario can occur provided that the 
neutron star fulfills some conditions. At first, the radius of the 
transition region has to lay inside the light cylinder radius of 
the neutron star, i.e Pa < Rlc = cP/2n, where P = lPi s is 
the rotational period of the neutron star, and c is the velocity of 
light. The above condition is fulfilled when, 

Pi > OmAB A ^M\l n . (5) 

Therefore, only relatively slowly rotating neutron stars can be 
considered. 

At second, the rotational velocity of the magnetosphere at 
R A has to be lower than the keplerian velocity of the accreting 
matter. The rotational velocity, 

y rot = 2ttR a /P * 2.5 x 10 9 B 4 ^ M^ n /P l cms" 1 . (6) 
is lower than the keplerian velocity, 

y k = (GM NS /P A ) 1/2 * 6.8 x lO 8 ^ 2 ^ 7 cm s" 1 . (7) 
for the NS with the periods, 

Pi > 3.7 B^Nql 11 . (8) 

This last condition separates the population of NS in the pro- 
peller phase (lower periods) from these ones in the accretor 
phase (larger periods). Note, that discussed below scenario for 
acceleration of electrons can be applied also to the propeller 
stage of accretion onto the NS. However, in the propeller phase 
the radiation field created by the matter accreting onto the NS 
is not uniquely defined. A part of the matter accreting onto NS 
can not be at present stage of knowledge linked to the basic 
parameters describing the model. We have to introduce addi- 
tional factor which describes the amount of the matter which 
falls onto the NS surface to the total amount of the matter pen- 
etrating the inner NS magnetosphere (i.e accreting and expelled 
from the vicinity of NS). 

The third condition relates the radius of the transition re- 
gion, Pa, to the capturing radius of the matter from the stellar 
wind. It is determined by the balance between the kinetic en- 
ergy of the wind with its gravitational energy around the NS. 
The capturing radius is described by, 

R c = 2GM NS /vl » 3.7 x 10 10 v 8 2 cm, (9) 

where v w = 10 8 vs cm s _1 is the velocity of surrounding matter 
measured in respect to the NS. The accretion from the stellar 
wind occurs when P c > Pa, which happens for, 

B l2 < 2.8 x 10 3 v„ 7/2 M?f . (10) 



For likely parameters of the neutron stars (classical and mil- 
lisecond pulsars), this limit is not restrictive. Note that v w cor- 
responds to the velocity of the stellar wind and/or the velocity 
of the neutron star on its orbit around the massive star. This last 
velocity can be estimated from, vns = 2nD/T = 5 x 10 7 Z)2/Pio 
cm s _1 , where D = IOOD2P0 cm is the radius of the orbit in 
units of 100 solar radii P , and T = lOPio days is the period of 
the binary system in units of 10 days. It is typically lower than 
the free fall velocity at the Pa- Therefore, we neglected it when 
deriving Eq.[3] 

Neutron stars with the periods within the range defined by 
Eq. 5 and Eq. 8 are in the propeller phase, those ones with the 
periods above that given by Eq. 8 are in the accretor phase, 
and those ones with the periods shorter than given by Eq. 5 
are in the ejector phase. As already noted, only in the accretor 
phase the accretion rate can be directly linked to the thermal X- 
ray emission from the NS surface. In the propeller phase, the 
amount of the matter which accrete onto the NS surface can 
not be uniquely determined from the basic parameters of the 
model. 

2.1. Acceleration of electrons 

In the conditions expected for the transition region (strongly 
magnetized and very turbulent medium), particles should be 
efficiently accelerated. In this paper we consider only acceler- 
ation of electrons. The energy gain rate of electrons with en- 
ergy E (and the Lorentz factor y) is often parametrized by the 
Larmor radius of electrons and so called acceleration parame- 
ter, 

P acc = {cE/r L * 2.6 x 10 4 f_!M^ 7 B- 2 5/7 erg s" 1 , (11) 

where £ = 10 _I £_i is the acceleration parameter, c is the ve- 
locity of light, = E/eB A is the Larmor radius, and e is the 
electron charge. The acceleration parameter contains all the un- 
known details of the acceleration process. Since the plasma 
moves in the transition region with the velocity which takes 
a significant part of the velocity of light, it seems proper to 
consider the values of ^ in the range ~ 0.1 - 0.01. During the 
acceleration process, electrons also suffer energy losses on the 
synchrotron process and the inverse Compton scattering of ra- 
diation from the massive star and the surface of the neutron star. 
These energy losses determine the maximum energies of ac- 
celerated electrons since their Larmor radius is typically much 
smaller than the characteristic dimensions of the considered 
scenario, i.e. Pl < Pa- This last condition allows in principle 
acceleration of electrons up to E e ~ 2M 4/ ) 7 P 12 8/7 PeV. However, 
as we show below, energy losses of electrons limit their maxi- 
mum energies significantly below this value. 

Electrons lose energy on the IC process in the Thomson 
(T) and the Klein-Nishina (KN) regimes. In general, this pro- 
cess can occur on the stellar radiation and on the thermal X-ray 
radiation from the surface of the polar cap on the NS. Let us 
estimate the photon energy densities from the star (p*) and the 
polar cap (p cap ) at the acceleration region, 

P* = ^t(tT ) 2 * 6 - lxl03? t(f f er 8 cm ~ 3 > ( 12 ) 
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Fig. 2. The maximum energies of accelerated electrons (a) and 
the maximum power which can be transfered to these electrons 
(b) as a function of the accretion rate of the matter from the 
massive star onto the neutron star. The results are shown for the 
neutron stars which have the surface magnetic field B — 10 12 G 
(solid line), and B = 10 9 G (dashed). The acceleration parame- 
ter is equal to £_i = 1. The estimated Fermi LAT sensitivity at 
1 GeV (1 yr) for the source at the distance of 1 kpc is L 7 ss 10 32 
ergs"' 



1 (thin dot-dashed line) and CTA sensitivity (50 hrs) at 20 
1 (thin dotted line). 



GeV Ly « 3 x 10 32 erg s 



(where r* = 3 x 10 4 r 4 K and cr is the Stefan-Boltzmann con- 
stant) and, 



Pcap 



1.4 x 10 8 M." /7 B" 8/7 erg cm" 3 , 



(13) 



where r cap = lO 7 ^ K is the temperature of the polar cap on 
the NS surface, R cap is the polar cap radius on the NS surface 
on which the matter accrete. In Eq. [13] we have assumed that 
the observed X-ray thermal emission is re-radiated from the 
region of the polar cap according to Lx = 7rR 2 ap crr 4 ap . In fact, 
the emission from the polar cap region is well described by the 
Black Body spectrum (see e.g. Zane et al. 2000). The radius 
of the polar cap region on the NS surface, on which the matter 
falls and from which the thermal X-ray emission is emitted, can 
be estimated from (assuming dipole structure of the magnetic 
field), 



^cap - (^ns/^a)' /2 



5 x 10 4 B7 2/7 M!f 



Then, the surface temperature of the polar cap is, 



r cap = (L x /^ 2 ap( r) 1/4 * 4.7 x 10 7 fl^ 7 M 5/28 



16 



K. 



(14) 



(15) 



Note that, the energy density of radiation from the polar cap 
at the distance of the turbulent region (Ra) can be explicitly 
expressed by the basic parameters of the model. 

Let us also estimate the energy density of the magnetic field 
at the transition region (Ra given by Eq.0), 



p B = B 2 J8n * 10 7 M 2// fi7? erg cm 



(16) 



The energy losses of electrons for every considered process 
(synchrotron and IC in the T regime) can be calculated from, 



Pi oss = (4/3)«r T py 2 * 2.7 x 10" ,4 p (cap , B ,*)y 2 erg s"\ (17) 
where <x T is the Thomson cross section. 



The energy losses of electrons on the radiation from the 
NS surface can dominate over energy losses on other targets 
only at low energies, i.e. when the scattering occurs in the 
T regime which happens for the Lorentz factors of electrons, 
Tt/kn < ^eC 2 /(3/tr*) * 20077 1 . In the KN regime, the en- 
ergy losses can be approximately estimated by introducing into 
Eq-Efrhe value of the Lorentz factor corresponding to the tran- 
sition between T and KN regimes, y^^y- For the range of pa- 
rameters defining the model, the IC energy losses of electrons 
in the KN regime becomes lower than the synchrotron energy 
losses for the Lorentz factors about ten times larger than T^/kn 
(typically at ~GeV energies). Therefore the maximum energies 
of accelerated electrons are determined by the balance between 
energy gains from the acceleration process (Eq.fTTTi and energy 
losses on synchrotron process (Eq. fTTT i, 



(18) 



r^-sxio 5 ^ 14 ^ 7 - 

It is clear that electrons can reach even TeV energies (see 
Fig. |2k) for realistic parameters of the model. However, TeV 
y-ray production will be strongly suppressed due to dominant 
synchrotron energy losses of electrons with TeV energies. 

We can also estimate the characteristic energies of syn- 
chrotron photons which can be expected in such a model by 
applying the approximate formula, 

e x * m e c 2 (B A /B cr )y 2 nax * L9£_, MeV, (19) 

where ra e is the electron rest mass, B CI = 4.4 x 10 13 G is 
the critical magnetic field, and the Lorentz factor of electrons 
ymax is given by Eq. [18] We conclude that synchrotron emis- 
sion can extend through the soft y-ray energy range, it should 
be detectable by satellites sensitive at these energies, e.g. the 
INTEGRAL observatory. 

Note that the bremsstrahlung energy losses of relativistic 
electrons in the matter inside the transition region can be ne- 
glected since their energy loss rate is relatively low, 

P br = m e c 3 p H y/X * 6.5 x 10" 5 pny erg s" 1 , (20) 

where X = 62 g cm" 2 is the radiation length for hydrogen, 
and pe = 10 17 pi7 cm" 3 is the density of matter. We can neglect 
the bremsstrahlung energy losses in respect to other energy loss 
rates since typical density of matter in the transition region is 
P 

cretion rate and the location of Ra). 



1.5 x 10 13 B 12 6/7 Mjg /7 cm" 3 (estimated for the known ac- 



2.2. Energetics 

The maximum power available for acceleration of electrons is 
limited by the energy extracted from the rotating neutron star 
by the in-falling matter. This matter from the stellar wind has 
to be accelerated to the velocity of the magnetic field lines at 
Ra in order to farther accrete onto the neutron star surface. The 
power which has to be transfered from the rotating NS to the 
matter can be estimated from 



M, 



acc^rot/^ 



3 X 10 M B?Mr Pf ergs" 1 . 



'12 '"16 'I "6 3 ■ ( 21 ) 

By using the limiting period given by Eq. [8] we get the upper 
limit on the available power in the accretor stage, 



L acc < 2.2 x 10 33 B 



33r-4/7 m 9/7 



12 



16 



erg s 



(22) 
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This power is shown in Fig. 2b for two example surface mag- 
netic field strengths of the neutron star (classical NS: B = 
3 x 10 12 G and millisecond NS: B = 10 9 G) as a function of 
the accretion rate. We assume that only a part, rj, of this power 
(applying typically rj = 0.1) can be converted to relativistic 
electrons in the turbulent and magnetized plasma at R\. Note, 
that for reasonable accretion rates this power is greater than the 
minimum required power of the GeV y-ray source which can 
be detected by the Fermi LAT telescope and possibly also by 
the planned next generation Cherenkov telescopes system CTA 
whose sensitivity should be about one order of magnitude bet- 
ter than presently available. 

The model discussed above for the accretor stage of the 
binary system can also operate for the propeller stage in which 
only a part of the matter can eventually accrete onto the surface 
of the NS. However, in this case it is difficult to estimate the 
power which can be transfered to relativistic particles due to 
the unknown observational signatures of the amount of matter 
which accumulates close to the shock region. 

In summary, electrons can be accelerated even to TeV ener- 
gies in the case of neutron stars inside X-ray binaries with rela- 
tively weak thermal X-ray emission from its surface in contrary 
to the powerful X-ray binaries in which only a few ten GeV 
electrons are expected. However, the situation is the opposite 
concerning the maximum power transfered to these relativistic 
electrons. Sources in which electrons are accelerated only to 
GeV energies should be powerful enough to be observable by 
the y-ray telescopes. 

3. Production of radiation 

We assume that electrons are accelerated in the turbulent transi- 
tion region with the power law spectrum to maximum energies 
estimated from Eq. [18] Eventually, the spectrum of electrons 
can be strongly peaked at the highest possible energies due to 
the synchrotron energy losses during acceleration process (the 
so called pile-up mechanism, see e.g. Protheroe 2004). As we 
have shown above, electrons lose energy on different radiation 
mechanisms but the dominant ones are: the synchrotron pro- 
cess (dominates at the highest energies); and the ICS of thermal 
radiation from the NS polar cap (contributes to the high ener- 
gies in the KN regime and dominates at lower energies in the T 
regime). We neglect the production of y-rays by electrons in the 
acceleration region in the scattering of stellar radiation since its 
energy density can be safely neglected in respect to the energy 
density of the magnetic field (at high energies) and the energy 
density of polar cap radiation (at low energies) (see Sect. 2.1). 
In order to check whether electrons can lose efficiently energy 
already in the transition region, we estimate their convection 
time scale with the matter falling onto the NS surface along the 
open magnetic field lines on, 

Tconv = ^A/Vrot = Pl/2n S. (23) 

Let us compare this time scale with the radiation time scale, 
which (for electrons with low enough energies) is determined 
by the IC losses (in T regime) in the radiation field from the 
polar cap of the NS, 

T mA ^m e ylP mi ^Q.2M- l l yn B%\- 1 s (24) 



From this comparison, we estimate lower limit on the Lorentz 
factor of electrons which lose efficiently energy before being 
convected onto the NS surface on, 

y=\3M-£ n B\%p- l l . (25) 

It is concluded that relativistic electrons lose their energy close 
to the acceleration place. 

In order to calculate the y-ray spectra produced by electrons 
inside the transition region, we simulate the energy loss process 
of electrons in the radiation field of the polar cap and at the 
magnetic field of the transition region by applying the Monte 
Carlo method. It is assumed that the magnetic field is very tur- 
bulent in the transition region (the distribution of electrons is 
isotropic) and that the radiation field comes from the region of 
the polar cap. For some range of model parameters, IC y-rays 
can be farther absorbed in the radiation field of the polar cap. In 
this way the IC e ± pair cascade is initiating. In order to deter- 
mine the conditions for the cascade process, in the next Section 
we calculate the optical depths for y-ray photons in the radia- 
tion of the polar cap. We calculate also the optical depths for 
y-rays in the radiation field of the massive companion in order 
to check whether this y-ray absorption process should be also 
taken into account when evaluating the y-ray spectra escaping 
from the binary system. 

4. Optical depths for y-rays 

As we have argued above, the cooling of electrons in such 
strong magnetic and radiation fields occurs very efficiently. 
High energy y-ray produced in the ICS of thermal photons from 
the polar cap can be also absorbed in this same radiation field 
and, in principle, also in the radiation field of the close mas- 
sive companion star. Below, we calculate the optical depths for 
the y-ray photons in these two radiation fields. We show that 
for specific conditions (determined mainly by the accretion rate 
onto the NS and the parameters of the NS), the y-ray spectrum 
emerging from the binary system is in fact formed in the cas- 
cade process occurring in the radiation field of the polar cap. 

4.1. Polar cap radiation 

We calculate the optical depths for y-rays assuming that they 
are produced by electrons in the transition region at the distance 
Ra from the stellar surface. In these calculations we assume the 
diluted black body spectrum for the polar cap emission with 
the dilution factor at the production site of y-rays estimated 
from (R C! l P /Ra) 2 - The average optical depths (averaged over the 
isotropic injection of y-rays) for different accretion rates onto 
the neutron star (i.e equivalent to different X-ray luminosities, 
see Eq. 1) are shown in Fig. [3] It is clear that for low accretion 
rates (< 10 16 g s _1 ), when electrons can be accelerated to ~TeV 
energies, the optical depths for y-ray photons are relatively low, 
i.e. they are below unity for energies above ~ 10 GeV. In this 
case GeV- TeV y-rays can escape from the radiation field of the 
polar cap without significant absorption. 

On the other hand, for large accretion rates (> 10 16 g s _1 ), 
electrons are accelerated at most to a few tens of GeV. y-rays 
produced by these electrons should be efficiently absorbed in 




log (E / GeV) 

Fig. 3. The optical depths for y-rays (averaged over the 
isotropic injection of y-rays) in the thermal radiation from the 
polar cap on the NS surface as a function of y-ray energies 
and for different accretion rates of matter onto the polar cap: 
M acc = 10 18 g s _1 (dotted curve), 10 17 g s~' (dashed), 10 16 g 
s _1 (solid), 10 15 g s -1 (dot-dashed), and 10 14 g s -1 (triple dot- 
dashed). The neutron star has the surface magnetic field equal 
to B NS = 3 x 10 12 G (upper figure) and 10 9 G (bottom). The 
distance of the acceleration region from the polar cap is 
defined by the above parameters. The thin dotted line marks the 
optical depth equal to unity. 

the radiation of the polar cap. Therefore, in this case the escap- 
ing y-ray spectra can be only obtained by calculating compli- 
cated IC e ± pair cascade in the radiation of the polar cap. We 
conclude that for large accretion rates y-rays, produced in this 
model, could be observable by the satellite telescopes {Fermi 
LAT and possibly AGILE). 

4.2. Massive star radiation 

The optical depths for y-ray photons injected at an arbitrary 
distance from the surface of the massive star have been calcu- 
lated for the first time in the general case (including also dimen- 
sions of the star) by Bednarek (1997, 2000). In fact, the optical 
depths for y-rays close to the surface of other stars can be eas- 
ily re-scaled from those early calculations. For example, y-ray 
photons with energies, E° propagating at specific distances D 



Fig. 4. The optical depths as a function of y-ray energies are 
shown for the case of the massive star in the massive bi- 
nary IGR J19140+098 discovered by INTEGRAL telescope: 
R* = 30R@, = 2.8 x 10 4 K and two distances from the 
star D = 3/?* (a) and 2/?+ (b). Specific curves show the optical 
depths for the injection angles a = 0° (thick solid curve), 30° 
(thick dashed), 60° (thick dotted), 90° (thick dot-dot-dashed), 
120° (thin solid), 150° (thin dashed), and 180° (thin dotted), 
measured from the direction towards the star. 

and at directions (defined by the angle a) close to the star with 
specific parameters (T a and R Q ) are related to the optical depths 
around arbitrary star with and in the following way, 

E° 

t(E* = -t,T*,R ir ,D,a) = S*S r t(E° T ,R ,D,a) (26) 
'At 

where Sj = T+/T , Sr = Ri,/R , distance D from the stel- 
lar surface of specific star is measured in the stellar radii. The 
example calculations of the optical depths for the case of one 
specific massive star (present inside the binary system IGR 
J19140+098) are shown in Fig [4] These optical depths be- 
come larger than unity for the y-rays with energies above ~20 
GeV. As we show latter, y-rays produced in the considered 
model have typically energy below ~10 GeV. Their absorption 
in the radiation field of the massive star can be safely neglected. 
Therefore, the absorption effects of y-rays should not introduce 
any modulation of the y-ray signal with the period of the binary 
in the case of accreting neutron stars. This is in contrast to the 
TeV y-ray binaries supposed to contain ejecting neutron stars 
(strong pulsar wind prevents the accretion of stellar wind), see 
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e.g. the model considered by Sierpowska & Bednarek (2005). 
The eventual modulation of the y-ray signal with the period of 
the binary can be only produced by the change of the accretion 
rate with the distance of the neutron star on its elliptic orbit 
around the massive star. 

5. Gamma-rays from the vicinity of accreting 
neutron star 

In the considered model, radiation is produced by electrons ac- 
celerated to maximum energies estimated by Eq.[l8]in the syn- 
chrotron and IC process. The synchrotron energy losses dom- 
inate at the largest electron energies (close to E m - dx ) and the 
IC process can only dominate at lower energies due to the KN 
cross section. However, as we have shown above, for large ac- 
cretion rates, y-rays produced by electrons can be farther ab- 
sorbed in the thermal radiation from the NS surface. We have 
shown in Fig. [3] that the average optical depths for y-ray pho- 
tons in the thermal radiation from the NS polar cap region in 
order to envisage for what parameters the cascading effects can 
become important. Therefore, in these cases we have to con- 
sider the production of the y-rays in the IC e ± pair cascade 
with additional synchrotron energy losses of primary electrons 
and secondary e ± pairs. On the other hand, for small accretion 
rates, the cascade does not develop, y-rays are produced in this 
case only as a result of cooling of primary electrons. In order to 
follow the process of y-ray production, we developed the nu- 
merical code which simulate the cooling process of electrons 
taking into account not only the y-ray production in the IC pro- 
cess but also the synchrotron energy losses of primary electrons 
and secondary cascade e ± pairs. Note that the acceleration of 
electrons occurs inside the inner pulsar magnetosphere where 
the magnetic field is relatively strong and depends on the dis- 
tance of the turbulent transition region from the NS surface. 
Since the synchrotron process dominates over the IC process at 
the high energy part of the injected electron spectrum, the y- 
rays are rarely produced with energies comparable to energies 
of primary electrons. Our IC e ± pair cascade code with syn- 
chrotron energy losses gives both, the synchrotron and the IC 
spectra expected in such a model. 

Two models for injection of relativistic electrons are con- 
sidered: 

1 . Electrons injected with the power law spectrum up to the 
maximum Lorentz factor y max (see Eq. ITSt . as expected 
in the stochastic acceleration mechanism. The power law 
spectrum of primary electrons is normalized to a part, 77, 
of the kinetic power which is transfered from the rotating 
neutron star to the matter in the transition region. 

2. Electrons injected with the Lorentz factors y m ax- Such in- 
jection spectrum can give good approximation in the case 
of acceleration process occurring with large radiative en- 
ergy losses (the pile up mechanism at the end of accelerated 
power law spectrum of electrons, see e.g. Protheroe 2004). 

These two injection models are likely to be the limiting 
cases of the real acceleration process of electrons in the turbu- 
lent region of the matter accreting onto the magnetized neutron 
star. 



We investigate two general scenarios describing the cases 
of pure accretor phase (all the matter arriving to the transition 
region falls onto the NS surface) and the intermediate accretor- 
propeller phase (only part of the matter arriving to the transition 
region accrete onto the NS surface and the rest of it is expelled 
from the vicinity of the NS in the propeller mechanism). 

5. 1 . The accretor stage 

The scenario for the accretion process onto neutron star in the 
accretor stage is much better defined since the thermal radiation 
field from the NS surface (depending on the dimension of the 
hot spot and its temperature) is uniquely determined by the pa- 
rameters describing the model. Therefore, the synchrotron and 
IC spectra produced by electrons depends on a relatively small 
number of free parameters, i.e. the accretion rate onto NS, the 
surface magnetic field of NS, and the acceleration rate The 
observed X-ray to y-ray power depends additionally on the 
parameter, rjD 2 , which combines the energy conversion effi- 
ciency from the transition region to relativistic electrons, rj, and 
the distance to the source D — 1D\ kpc. At first, we performed 
calculations for the power law spectrum of injected electrons 
N(e) = AE~ S between E, mx and E m [ n , where A is the normal- 
ization coefficient equal to A — rjL acc {2 - s)/{E < ^ > - E®7^). 
Note, that A only weakly depends on E m \ n for the spectral in- 
dices not far from 2. In our example calculations, we apply the 
spectral index equal to 2.2 and choose E m \ n = 30 MeV since we 
cool electrons only to such minimum energies. Electrons with 
such energies are not able to produce y-rays detectable by the 
Fermi LAT. The calculated spectra are also shown for different 
values of M, £ and B (see Fig. [5j. 

In Figs. Hh-d, we consider the case of classical relatively 
young NS (B » 10 9 G), rotating with the limiting period de- 
fined by Eq. [8] When confronting with the sensitivity of the 
Fermi LAT telescope, it becomes clear that y-ray emission from 
such objects can be detected at energies below ~ 1 GeV in the 
case of large accretion rates (significantly above 10 16 g s _1 ). 
Moreover, the power law synchrotron spectra, corresponding 
to these accretion rates, should be also easily detected by the 
soft and hard X-ray detectors presently on the orbit. Note that 
due to strong magnetic field and large energies of primary elec- 
trons, the synchrotron spectra can extend in some cases also 
through the soft y-ray energy range. In Fig.[3£-f, we show the 
synchrotron and IC spectra for the case of the NS with low sur- 
face magnetic field, i.e. characteristic for so called millisecond 
pulsars. Also in the case of these slowly magnetized neutron 
stars, y-ray emission can be detected by the Fermi LAT tele- 
scope, provided that the accretion rate onto the NS is above a 
few 10 14 g s _1 . We conclude that y-ray telescopes presently on 
the orbit should be able to detect y-ray emission at a few hun- 
dred MeV from the class of accreting NS. These NS are char- 
acterized by the thermal emission from its surface at a level 
above ~ 2 x 10 36 erg s _1 in the case of strongly magnetized 
neutron stars, and above a few 10 34 erg s" 1 in the case of NS 
with the magnetic field characteristic for millisecond pulsars 
provided that the source is at the distance of 1 kpc. These limit- 
ing X-ray luminosities have been derived by using the relation 
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between the accretion rate and the thermal emission from the 
NS surface (see Eq.Q~|l. 

We also investigate how the detectability of such accreting 
NSs depends on the spectral index of the injected electrons. 
In Fig. [6] the synchrotron and IC y-ray spectra are shown for 
the range of spectral indexes 5=1-4 and for fixed other 
parameters describing the model: M - 10 17 g s _1 , £ = 0.1, and 
B — 10 12 G. Electron spectrum has to have spectral index lower 
than ~ 3 in order to produce IC y-rays above the sensitivity of 
Fermi LAT telescope. Such spectral indexes are expected in the 
stochastic models for particle acceleration in turbulent media. 

As already noted above, the acceleration of particles under 
strong radiative energy losses may result in their accumulation 
at the highest energies determined by the balance between en- 
ergy gains from the acceleration mechanism and energy losses 
(see e.g. Protheroe 2004). In such cases, most of the energy 
is accumulated in particles with the highest possible energies. 
At first approximation, we can consider that the spectrum of 
particles is mono-energetic. Therefore, in Fig.|7]we show also 
the photon spectra expected in the model for the case of mono- 
energetic injection of electrons. In such a case, IC y-ray spec- 
tra are flat. Their detectability is more difficult than in the case 
of injection of electrons with the power law spectra. This is 
due to the fact that electrons with maximum possible energies 
lose most of their energy on synchrotron process. Therefore the 
IC spectra are reduced and synchrotron spectra are enhanced 
(compare corresponding cases shown in Figs.|5]and|7]l. We con- 
clude that accreting NS, which have the highest chance to be 
detected by the Fermi LAT telescope, should have the interme- 
diate values of £ and B and large accretion rates (see Fig.[7J>d). 
Such parameters allow acceleration of electrons to energies of 
the order of a few tens GeV, i.e close to the energy region where 
ICS process starts to become comparably efficient to the syn- 
chrotron process. 

The y-ray photons produced by electrons in the transition 
region of the accretion flow onto the NS originate close to the 
massive companion star which creates strong radiation field. In 
Sect. 4.2, we show the optical depths for y-rays in the radiation 
field of the massive star with the parameters typical for mas- 
sive binary systems (the example case of the massive binary 
IGR J19140+0951). As we have shown above, primary elec- 
trons accelerated to the maximum energies cool at first mainly 
on the synchrotron process. Therefore y-rays produced by them 
in the IC process have energies typically below ~ 10 GeV. On 
the other hand, the optical depths in the massive star radiation 
(see Fig. HJ are quite large but at energies which are clearly 
above energies of y-rays produced in our model. Therefore, we 
conclude that in the case of considered here model the absorp- 
tion of y-rays in the stellar radiation can be neglected. 

In Fig. 7a,c,e, we also show the spectra of thermal radia- 
tion from the NS polar cap region for the example accretion 
rate M acc = 10 16 g s _1 and different values of the magnetic 
field on the NS surface (B = 10 12 G (a), 10 11 G (c), and 10 9 
G (e). These spectra clearly dominate over the nonthermal syn- 
chrotron spectra produced by relativistic electrons below a few 
keV. In fact, as we already noted above, such soft X-ray ther- 
mal excesses have been reported from some INTEGRAL hard 
X-ray massive binaries. In reality this soft thermal X-ray emis- 



sion may be significantly modified due to the interaction with 
the matter accreting onto the NS surface, the matter of the stel- 
lar wind, and the matter onto the surface of the close massive 
companion. In our calculations we are not able to take these 
possible modifications into account. 

5.2. The intermediate accretor-propeller stage 

As noted above, in the pure propeller stage the matter arriving 
to the transition region at the Alfven radius, R\, is expelled by 
the centrifugal source from the vicinity of the neutron star pref- 
erentially along the rotational axis. Here we consider the inter- 
mediate case in which a small part of the accreting matter is 
able to penetrate onto the NS surface but most of it is expelled 
outside the neutron star. It is not clear whether such process can 
last stationary in time or it occurs only non-stationary when the 
accretion rate onto the neutron star changes in time e.g. due to 
the clumpy wind or the elliptic orbit of the NS around the mas- 
sive star. In order to take into account the effects of only partial 
accretion of matter onto the NS surface, we introduce the pa- 
rameter, k, which is the ratio of the matter accreting onto the NS 
surface to the whole amount of matter arriving to the transition 
region (i.e. the matter accreted onto the surface and expelled 
from the vicinity of NS). In such an intermediate case, the tem- 
perature of the polar cap region on the NS surface (given by 
Eq.l23l should be a factor k 1 ^ lower than estimated in the case 
of complete accretion. However, in principle the power avail- 
able for acceleration of electrons can be larger than expected 
from period limiting the accretor and propeller stage (given by 
Eq. [8]). Instead, it is rather limited by Eq. [3] In conclusion, as 
a result of only partial accretion of matter from the transition 
region onto the NS surface, the surface temperature of the polar 
cap is lower and the power transfered to electrons can be larger 
than expected for the pure accretor stage. 

In Fig. [8] we show the example y-ray spectra obtained in 
the case of the partial accretion of matter onto the NS surface 
assuming different accretion rates defined by the factor k. We 
also scale the power transfered to relativistic electrons by the 
factor k since more energy can be extracted from the transi- 
tion region in this case. Based on these calculations, we expect 
that the y-ray spectra produced in the intermediate stage are 
steeper, although they have larger luminosities. Therefore, they 
should be easier detected by the y-ray telescopes (see Fig.[8]for 
the comparison with the sensitivity of Fermi LAT telescope for 
the source at the distance of 1 kpc). 

6. Discussion and Conclusions 

We propose that accreting neutron stars inside compact mas- 
sive binary systems produce y-ray fluxes which can be observ- 
able by the satellite telescope Fermi LAT. To show this we con- 
sider acceleration of electrons at the turbulent, strongly mag- 
netized, transition region in the inner neutron star magneto- 
sphere which appears as a result of the interaction of accreting 
matter with rotating NS magnetosphere. Relativistic electrons 
produce X-rays and y-rays as a result of synchrotron and IC 
processes occurring on the thermal radiation from the NS sur- 
face. The cooling process of electrons is followed by applying 
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the Monte Carlo method since some of produced y-rays can be 
farther absorbed in the thermal radiation. We showed that syn- 
chrotron emission from primary electrons and secondary cas- 
cade e ± pairs can extend up to the MeV energy range and the 
IC emission can appear up to a few GeV. 

A physical realization of such a scenario can be charac- 
teristic for the newly discovered by the INTEGRAL observa- 
tory a class of obscured massive binary systems which show: 
(1) a hard, power law X-ray emission and, (2) evidences of 
soft X-ray black body component. As an example, we consider 
here in a more detail the high mass X-ray binary system, IGR 
J19140+0951, discovered by the INTEGRAL (Hannikainen et 
al. 2003). A compact object in this system is likely to be a neu- 
tron star (Cabanac et al. 2005). IGR J19140+0951 is at the dis- 
tance of ~ 2-3 kpc (Rahoui et al. 2008). The hard X-ray 
emission up to 100 keV with the spectral index 2.39 + 0.11 is 
observed from this source (Hannikainen et al. 2004). Its X-ray 
luminosity in the high state was ~ 3.7 x 10 37 x (D/lOkpc) 2 erg 
s _1 (Rodriquez et al. 2005). Unfortunately, the basic param- 
eters of the neutron star in IGR J 19 140+0951 are unknown. 
Therefore, we investigate the range of model parameters which 
are consistent with the hard power law spectrum in the X-ray 
energy range by the INTEGRAL. Comparison of the observa- 
tions with the example calculations are shown in Fig. [9] It is 
clear that for some parameters the y-ray flux predicted at a few 
hundred MeV can be detectable by the extensive Fermi LAT 
observations. 

We performed calculations for the neutron stars with the 
parameters characteristic for the "classical" radio pulsars (sur- 
face magnetic field strength B ~ 10 12 G and periods of the 
order of seconds) and for the millisecond pulsars (B ~ 10 9 G, 
periods of a few to several milliseconds). It is clear that in or- 
der to produce y-ray fluxes observable by the Fermi LAT tele- 
scope, the neutron stars in binary systems at the distance of a 
few kpc should collect the matter from the wind at a relatively 
large accretion rates (> 10 16 g s ). However, in the case of 
millisecond pulsars these accretion rates can be significantly 
lower (above a few 10 14 g s _1 ). Therefore, in principle, the mil- 
lisecond pulsars seems to be more favorite y-ray sources. On 
the other hand, millisecond pulsars have only low mass com- 
panions, so the large accretion rates are not expected in such 
binary systems. 

Note, that synchrotron X-ray and IC y-ray emission is pro- 
duced in such a model almost isotropically. As we have shown 
above, we do not expect any modulation of the y-ray signal 
with the period of the binary system due to the selective ab- 
sorption in the companion star soft radiation, as it is expected in 
the case of TeV y-ray binary systems (e.g. LS 5039 and LSI 61 
303). However, the modulation might be related to the change 
of the accretion rate in the case of elliptic orbit of the neutron 
star. It is expected that the accretion rate should increase when 
the neutron star is closer to the companion star. Therefore, we 
predict that the largest fluxes of X-rays and y-rays should be 
observed closer to the periastron passage of the neutron star. 
Moreover, the X-ray and y-ray emission should be correlated 
(see spectra in Figs. 5-8). 

In the model considered here we do not take into account 
the scattering of nonthermal synchrotron radiation by the ac- 



celerated electrons. As we have shown in Fig. 7, the thermal 
radiation clearly dominates at energies below a few keV over 
the nonthermal synchrotron radiation produced by accelerated 
electrons and secondary cascade e ± pairs. Since their energy 
densities at the acceleration site scales in a similar way with 
the distance from the NS, the energy losses on thermal radia- 
tion has to dominate over the energy losses on nonthermal radi- 
ation. Therefore, we can safely neglect the cooling of electrons 
on this nonthermal radiation. Note moreover, that the scattering 
of nonthermal X-rays with energies above a few keV becomes 
inefficient since it occurs in the Klein-Nishina regime. 

Similar processes to these considered in this paper are also 
expected in the case of non-spherical accretion of matter onto 
the rotating neutron star. For example, in the case of disk ac- 
cretion, the corresponding accretion rate should be scaled by a 
part of sphere which is intercepted by the accretion disk. Then, 
the turbulent region on the inner disk age, in which accelera- 
tion of electrons occurs, is also limited to corresponding part 
of the sphere. Only the IC process may become more compli- 
cated in this case since thermal radiation from the hot spot on 
the neutron star surface, visible by relativistic electrons, can be 
partially obscured by the neutron star surface. 

In conclusion, we predict that present satellite telescopes 
can discover a new class of y-ray sources, i.e. accreting neutron 
stars inside binary system. The emission from these sources 
should be characterized by a strong hard, synchrotron spectra 
and IC y-rays limited to energies below ~ 1 GeV, due to the 
very strong synchrotron energy losses of relativistic electrons 
accelerated in the vicinity of neutron stars. 
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Fig. 5. The IC y-ray spectra (on the right) and synchrotron spectra (on the left) produced in a sequence of radiation processes by 
relativistic electrons. The IC process occurs on the thermal radiation field produced by the matter accreting onto the neutron star 
surface. The X-ray and y-ray spectra are shown for different accretion rates: M = 10 18 g s _1 (dashed curve), 10 17 g s _1 (solid), 
10 16 g s (dotted), M acc = 10 15 g s _1 (dot-dashed), 10 14 g s _1 (triple dot-dashed). The spectra are shown for the acceleration 
parameter £ = 0.1, the surface magnetic field of the NS Bns — 10 12 G (upper figures), £ = 0.01 and Z?ns = 10 n G (middle), and 
£ = 0.1 and Z?ns = 10 9 G (bottom). The spectrum of accelerated electrons is of the power law type {dN(E)/dE oc E~ s ) with the 
spectral index s = 2.2 between 30 MeV and the maximum energy given by Eq.[l8] The factor determining the power transfered 
from the transition region to relativistic electrons depends on the energy conversion efficiency, rj, and the distance to the source, 
D = IDi kpc. It is assumed to be equal to r/D^ 2 =0.1. The period of the pulsar is equal to the limiting value given by Eq.[8] The 
level of sensitivity of the the EGRET and GLAST telescopes are marked by thin dotted and double-dashed lines. 



12 W. Bednarek: Gamma-rays from accreting neutron stars 




'—6 -5 —4- -3 -2 -1 ' —1 O 1 

log (E / GeV) log (E / GeV) 



Fig. 6. As in Fig.Qbut for different spectral indexes of injected electrons: s — 4 (dashed), 3 (dotted), 2.2 (solid), 1 (dot-dashed). 
The parameters describing the model are the following: £ = 0.1, the surface magnetic field of the NS B^s = 10 12 G, and the 
accretion rate M acc = 10 17 g s~', . 
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Fig. 7. As in Fig. [5] but for mono-energetic injection of electrons with energies given by Eq. [18] The specific curves show the 
results obtained for these same parameters as in Fig. [5] The spectrum of thermal radiation emitted from the neutron star surface 
is also shown for the typical accretion rate M acc = 10 16 g s _1 and corresponding surface magnetic field strengths (thin dotted 
curves in (a), (c), and (e). 
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Fig. 8. As in Fig. [5] but for different rates describing the amount of the matter accreting onto the NS surface in respect to the 
whole amount of the matter arriving to the transition region (accreting and expelled from the vicinity of the NS in propeller 
mechanism): k — 1 (solid curve), 0.1 (dashed), and 0.01 (dotted). The basic parameters of the model are the following: M = 10 17 
gs" 1 ,B= 10 12 G, and£ = 0.1. 
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Fig. 9. (a) The comparison of the hard X-ray spectra detected by the INTEGRAL observatory from the massive binary IGR 
J 19 140+0951 with the synchrotron spectra expected from the vicinity of neutron star for the pure accretor model with the 
parameters: M = 2 x 10 17 g s _1 , B - 10" G, and £ = 0.1, spectral index of electrons s = 2.8 (solid curve), M = 3 x 10 17 g s _1 , 
B = 10 9 G, and £ = 0.1, s = 4 (dot-dashed), the intermediate accretor-propeller models for: M = 10 17 g s _1 , B = 10 11 G, £ = 0.1, 
and k = 0.1 and s = 2.8 (dotted), and M = 10 17 g s"\ B = 10 12 G, £ = 0.1, and k = 0.01 and s = 2.8 (dashed), (b) The IC 
y-ray spectra expected in the model for the parameters mentioned above based on the normalization of the synchrotron spectra 
to the INTEGRAL observations in the X-ray energy range. The level of sensitivity of the the Fermi LAT telescope is marked by 
a double-dashed line. 



